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Abstract 

We study the isospin singlet pairing in asymmetric nuclear matter with 
nonzero total momentum of the condensate Cooper pairs. The quasiparti- 
cle excitation spectrum is fourfold split compared to the usual BCS spectrum 
of the symmetric, homogeneous matter. A twofold splitting of the spectrum 
into separate branches is due to the finite momentum of the condensate, the 
isospin asymmetry, or the finite quasiparticle lifetime. The coupling of the 
isospin singlet and triplet paired states leads to further twofold splitting of 
each of these branches. We solve the gap equation numerically in the isospin 
singlet channel in the case where the pairing in the isospin triplet channel 
is neglected and find nontrivial solutions with finite total momentum of the 
pairs. The corresponding phase assumes a periodic spatial structure which 
carries a isospin density wave at constant total number of particles. The 
phase transition from the BCS to the inhomogeneous superconducting phase 
is found to be first order and occurs when the density asymmetry is increased 
above 0.25. The transition from the inhomogeneous superconducting to the 
unpaired normal state is second order. The maximal values of the critical 
total momentum (in units of the Fermi momentum) and the critical density 
asymmetry at which condensate disappears are Pc/pf = 0.3 and Uc = 0.41. 
The possible spatial forms of the ground state of the inhomogeneous super- 
conducting phase are briefly discussed. 
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I. INTRODUCTION 



The theory of fermion pairing when the fermions, which are bound in Cooper pairs, 
he on different Fermi surfaces was addressed in the early 1960s in the context of metalhc 
superconducting alloys containing paramagnetic impurities Recently, the theory has 

received renewed attention in several contexts, including the pairing in asymmetric nuclear 
matter [^-11] and color superconductivity in flavor-asymmetric high-density QCD |12,13 • 



This paper elaborates on an earlier suggestion ^ that the ground state of the supercon- 
ducting asymmetric nuclear matter at large asymmetries corresponds to a pair condensate 
with nonzero total momentum of the Cooper pairs. The argument is based on the observation 
by Larkin and Ovchinnikov JI^ and Fulde and Ferrell who first showed, in the context 
of the metallic superconductors, that the Bardeen-Cooper-Schrieffer (BCS) equations admit 
solutions with nonzero total momentum of Cooper pairs. In the configuration space such 
a condensate forms a periodic lattice with finite shear modulus. The resulting spatially in- 
homogeneous superconducting state is called the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) 
phase. 

The occurrence of pair correlations crucially depends upon the overlap between the neu- 
tron and proton Fermi surfaces; the pairing gap is largest in the isospin symmetric case and 
is suppressed as the system is driven out of the symmetric state. The thermal smearing 
of the Fermi surfaces promotes the pairing, but, however is ineffective when the separation 
between the surfaces is large compared to the temperature. If the total momentum of the 
Cooper pairs is zero, the Fermi surfaces (for homogeneous systems) are located on concen- 
tric spheres. If, however, a Cooper pair moves with a finite momentum, the centers of the 
spheres are shifted. This allows for the pairing, even for vastly different radii of the spheres, 
since now the nonconcentric spheres may intersect. The overlap regions then provide the 
kinematical phase space for pairing phenomena to occur. 

The microscopic calculations, based on the BCS theory for the bulk nuclear matter show 
that the isospin-symmetric matter supports Cooper-type pair correlations in the ^Si-^Di 
partial-wave channel due to the tensor component of the nuclear force. The energy gap 
is of the order of 10 MeV at the empirical saturation density (Refs. [§-[n| and references 
therein) if one assumes that the effective pairing interaction can be approximated by the bare 
interaction and if the renormalization of the single particle spectrum reduces the density of 
the state at the Fermi surface by a factor of 2. 



There is little evidence for large gap isospin singlet pairing in ordinary nuclei []18| , which 
is evidently suppressed due to the spin-orbit splitting The laboratory data do not 

exclude the possibility that the bulk nuclear matter, as encountered in the supernovas and 
neutron stars, may support large gap pairing in the isospin singlet channel. In the model of 
"nucleon stars" the kaon condensation implies nearly isospin symmetric matter, in which 
case the isospin singlet pairing can play a major role in determining the cooling and rotation 
dynamics of such objects. Note, however, that in the models without meson condensates 
the proton concentration in supernova and neutron star matter is of the order of 5% — 30% 
and these asymmetries are too large to allow for neutron-proton pairing. In the high- 
density regime, the hyperon-rich neutron star matter may be much more symmetric than at 
the densities around the saturation density and therefore can support neutron-proton 
pairing, most likely due to the attractive partial- wave interaction M. If the nucleon- 
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hyperon and hyperon-hyperon interactions are attractive, the pairing among fermions lying 
on different Fermi surfaces, and in particular the formation of the LOFF phase, could be a 
realistic possibility in hyperon rich matter. Finally, in superstrong magnetic fields of highly 
magnetized neutrons stars (magnetars) the Pauli paramagnetism will cause a splitting in 
the Fermi energies of the spin-up and spin-down fermions; in this case the pairing in the 
isospin triplet channel among the S-wave paired neutrons should exhibit the properties of 
the LOFF phase.Q 

Before turning to the main body of the paper, we draw the reader's attention to Ref. |[T6[ 



who considered the color superconducting LOFF phase in the context of high-density QCD 
and Ref. who considered the finite momentum pairing between quarks in the particle- 
hole channel (chiral condensate). Our work, to some extent, parallels the former reference, 
but we do not attempt any comparison at this stage, as the formalisms and contexts are 
entirely different. 

In Sec. II we derive the BCS equations, which include the effects of the finite momentum 
of the Cooper pairs, within the finite temperature real-time Green's functions formalism. 
The numerical solutions of these equations are shown in Sec. III. Sec. IV contains a summary 
of the results. 



II. FORMALISM 

Below, we shall use the real-time Green's functions extended to the Nambu-Gor'kov 
space to account for pair correlations. The single particle retarded Green's function in this 
space is defined as usual: 

"'^'"'''''^-'[fUx^x^) Gi,{x„X2))-[(T,iji{x,)4{x2)) {t,4{x,)Mx2)) ) ' 

where a and (3 stand for discrete quantum numbers (spin, isospin, etc.), ip"^ and are the 
nucleon creation and annihilation operators, x = {r,t) denotes the space-time coordinate, 
and Tf is the time-ordering symbol. The averaging is carried out over the equilibrium 
ensemble at a fixed density and temperature. The equation of motion for the matrix Green's 
function is given by the time-dependent Dyson equation 

G~^{Xi)Gap{Xi,X2) = iSapS{Xi - X2) + « ^ / rf^Xg S^^ (Xi , Xg)^^/? (Xg , X2) , (2) 

7 

where 1 is a unit matrix in the Nambu-Gor'kov space, the inverse free-particle propagator 
and the self-energy matrixes are 

f^-i(^_(G-\x) \ y / S,^(xi,X2) A„^(xi,X2) \ 

^"^-[ [G-(x)]* J ' ^"/^(^-^^) = Ul,(x„X2) ^Ux„X2) ) ' 



This problem is currently under study. 
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where G~^{x) = id/dt + V^/2ma + fia, fia is the chemical potential, and ma is the bare 
mass. The self-energy matrix is defined according to the rules of the usual diagram technique 
in the terms of G and the four-fermion interaction vertex T. In particular, the anomalous 
self-energy, which incorporates the pair correlations, is given by 

In the following we shall be interested in stationary (time-independent) and spatially in- 
homogeneous solutions of the equations above in the quasiclassical approximation. This 
approximation holds when the characteristic length scales of the spatial variations of the 
macroscopic condensate are much larger than the inverse of the momenta involved in the 
problem ~ Pf, where is the Fermi momentum. The quasiclassical counterparts of the 
equations above are obtained by going over to the center of mass X = {xi + X2)/2 and 
relative x = Xi — X2 coordinates in the two-point functions and carrying a Fourier trans- 
form with respect to the relative coordinates: G{x,X) G{u,p, R,T), where u,p are the 
relative frequency and momentum, and X = {R,T). As the variations of the propaga- 
tors and self-energies are slow on the scales of the order of R, keeping the leading order 
terms in the gradient expansion is accurate to order ~ 0[{pfR)~^]- Carrying out a Fourier 
transformation with respect to R, we arrive at the Dyson equation for the quasiclassical 
functions: 

where 



e^^ = (P/2 ± pf jlrria - /i^ ± Re S^/? - Im S^/?, (6) 

and Re S^^ = (S^/j — Sjj^)/2, Im S^/? = (S^/? + S|j,^)/2; all propagators and self-energies 
are functions of lo^p and P (the dependence on center-of-mass time is dropped in the 
stationary limit). Equation (1) is a (4 x 4) matrix in the spin-isospin space in general. The 
number of degrees of freedom can be reduced since the fermionic wave function of paired 
fermions must be antisymmetric. In the case of spin and isospin conserving forces the normal 
Green's functions and self-energies are diagonal in the spin and isospin spaces. It is sufficient 
to consider the anomalous propagators, e.g., in the isospin space, since the resulting spin 
structure, for S wave interactions, is uniquely determined for each isospin combination. The 
quasiparticle excitation spectrum is determined in the standard fashion by finding the poles 
of the propagators in Eq. (^: 



^±± = ± + ^Tr(AAt) ± i^[Tr(AAt)]2 - 4Det (AAt). (7) 

Here A = Aq,/}, es = (e+ + e^)/2, and eA = (e^~e^)/2. The new quasiparticle spectrum has 
four branches. The fourfold splitting of the BCS spectrum is due to (a) isospin asymmetry 
and/or the finite momentum of the condensate and/or the finite lifetime of the quasiparticles; 
and (b) the coupling of the pairing gaps in different isospin channels. If we restrict ourselves 
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to the neutron-proton pairing in the ^Si-^Di channel, which is justified when A„„, App <ti 
Anp, then Aa/s = cTxA {ax is the first component of the vector of Pauh matrixes). The 
spectrum, in this case, simphfies to 



u;± = e^±v/e| + |A|2, (8) 

where the symmetric and asymmetric parts of the spectrum (which are even and odd with 
respect to the time-reversal symmetry) are defined as 

es = P^/8m + p^/2m + Re S - /i, eA = P ■ p/2m + Im S - (9) 

Here /i = (/i„ + /ip)/2, 5fx = (/i„ - fip)/2 and Re S = (S„„ - Stp)/2, Im S = (S„„ + Sjp)/2 
(subscripts n and p refer to protons and neutrons). The limit ^ corresponds to the I3CS 
pairing in the isospin symmetric nuclear matter. It is explicit now that the spectrum (H) 
is twofold split due to three factors, the isospin asymmetry {6fi ^ 0), the finite- momentum 
of the Cooper pair {P ^ 0), and the finite lifetime of the quasiparticle (ImS ^ 0). Below, 
we shall keep the first two factors, and shall neglect the last one, since the quasiparticle 
approximation is valid in the density-temperature range of interest (densities around the 
nuclear saturation density and temperatures ~ 10 MeV). The solution of the Dyson equation 
(B) is 



Gn/p = + (10) 




F = UpVp — , (11) 

\UJ — UJ+ + IT] UJ — UJ- + IT] ' 



where the Bogolyubov amplitudes are 



Let us turn to the solution of the gap equation (§), which in the quasiclassical limit takes 
the form 

A(p,P) = 2 j ^0^Tico',p,p',P)lmF{co',p',P)fico'), (13) 

where /(a;) = [exp(/3a;) -|- l]~^is the Fermi distribution function and /5 is the inverse tem- 
perature; the effective pairing interaction F is assumed real and will be replaced by the bare 
interaction below. The u integration is straightforward in the quasiparticle approximation, 
since the frequency dependence of the propagator is constrained by the on-shell condition. 
Further progress requires partial wave decomposition of the interaction, which can be done 
after an angle averaging in the remainder functions on the right-hand side of the Eq. (|13|). 
The result of this procedure is 

Az(p,P) = -E / ^^/'(P.P')^=#P— ([/(^+)-/(^-)]), (14) 
(2^) 2Je| + A(p',P)2 
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where (...) denotes the average over the angle between the relative and total momenta, and 
A{p, Py = Aq{p, PY + A2(j9, PY is the angle-averaged gap. Here the pairing interaction is 
approximated by the bare neutron-proton interaction V{p,p') in the '^S'l-'^Di-channel. The 
self-consistent procedure of the determination of the gap function requires a normalization 
condition for the net density p = pn + Pp of the system at a fixed temperature and the 
magnitude of the total momentum P. The corresponding expression is provided by 

P^/,{P) = -2Y: J ^ImG„/,(^,p,P)/M = E/ (0jH/(^±) + ^'/(^t)}, (15) 

where a stands for quasiparticle spin and the second equality follows in the quasiparticle 
approximation. The coupled equations (0) and ([TsD should be solved simultaneously. 

To find the true ground state we have to minimize the free energy of the system at fixed 
total density and temperature. In the mean-field approximation the entropy of the system 
is given by the combinatorical expression 

S = -2kB E{/(^+) ln/(^+) + />+) ln/>+) + /(^-) ln/(^_) + />_) In />_)}, (16) 

where f{uj±) = [1 — f{uJ±)], and ks is the Boltzmann constant. The internal energy, defined 
as the grand canonical statistical average of the Hamiltonian, U = {H — fi^^'^Pn — p^'^'^pp), 
reads 



2 / [^'■^n(p) + e-n,{p)] + E VwiP^P) '^i{p)Mp)], (17) 



nn/p{p) = ulf{uj±) + fp/(^^), iy{p) = UpVp{[f{uj+) - f{uj-)]). (18) 



U 

where 



The first term in Eq. ([T7| ) includes the kinetic and non-pairing energies of the quasiparti- 
cles. The second term includes the BCS mean-field interaction among the particles in the 
condensate. The free energy (at fixed density and temperature) is defined as 

(F),,;3 = f/-r'5. (19) 

The true ground state of the system minimizes the free-energy difference {5F)p^p between 
the superconducting and normal states [the free energy in the normal state follows from Eqs. 
(0) and (|T3) when A = 0]. 



III. RESULTS 

The main focus of the numerical calculations shown below is the effects of the finite 
momentum of the Cooper pairs and the emergence of the LOFF phase in the asymmetric 
nuclear matter. A number of simplifying assumptions went into these calculations: first, the 
pairing interaction is approximated by the bare interaction; i.e., the effects of the screening 
of the pairing interaction are ignored. Second, we employ the quasiparticle approximation 
and set the effective mass of the quasiparticles equal to their bare mass. Third, we ignore 
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the coupling between the pairing in the isospin triplet and singlet states. The first two 
approximations change the absolute magnitude of the paring gap by affecting, respectively, 
the strength of the interaction and the density of states at the Fermi surface. To remove the 
dependence on the absolute scale of the gap we present the results normalized to the pairing 
gap in the symmetric matter at zero total momentum of the pairs. The third approximation 
is valid whenever the pairing in the isospin singlet channel is much larger than in the isospin 
triplet channel. This could be the case since the strength of the interaction in the ^Si-^Di 
in the free space is much larger the one in the ^5*0 channel and these channels are attractive 
in the same range of the energies. This argument, however, implicitly assumes that effects 
of the quasiparticle renormalization and the screening of the pairing interaction are of the 
same order in both channels, which could be false. 

Figure 1 shows the pairing gap A{pp) in the ^Si-^Di partial wave channel as a function 
of the isospin asymmetry, defined as a = (p„ — Pp)/p and total momentum P in units 
of the Fermi momentum. The pairing interaction has been approximated by a separable 
form of the Paris nucleon-nucleon interaction (PESTl of Ref. ||22|). The pairing gap is 



normalized to its value in the symmetric and zero-total- momentum case Aqo = 14.38 MeV 
at fixed total density p = 0.16 fm~^ and the temperature j3~^ = 3 MeV. The results reported 
here are relevant for the low-temperature regime (/3~^/Aoo = 0.208 <^ 1); the temperature 
dependence of the LOFF phase, in particular the T —>■ limit, will be discussed elsewhere. 

The absolute magnitude of the gap is consistent with the previous results based on the 
free-single-particle spectrum (although, note that the gap in Ref. is by 15% smaller, 
since there a rank 4 potential has been used instead of the rank 1 potential used in this 
work). A renormalization of the single particle spectrum, for example within the Brueckner 
theory, leads to a decrease of the gap by a factor of 2, see Ref. 0. This reduction also 
affects the critical asymmetry at which the BCS state disappears, by reducing it, e.g., at 
nuclear saturation density, from 0.35 for the free-particle spectrum to 0.11 for the Brueckner- 
renormalized particle spectrum . Therefore, the absolute magnitude of the asymmetry, 
at which the transition from the BCS to the LOFF phase occurs, and its critical value, 
at which the LOFF phase disappears, will be reduced by roughly a factor of 3, if the 
renormalization of the single particle spectrum is carried out within the Brueckner theory. 
For a = the gap is maximal at P = 0, decreases as the total momentum is increased, and 
vanishes at the critical total momentum Pc,o = 0.558pF. For P = the gap again decreases 
as a function of a and vanishes at ac,o = 0.37. The onset of the LOFF phase is signaled 
by the change of the shape of the constant a slices in the a-P plane: for a > 0.25 the 
maximum of the g function of P shifts from the P = to P 7^ values; i.e., the 

condensation energy becomes maximal for P 7^ 0. The maximum is located at P ~ 0.3pF 
and is not sensitive to the value of a > 0.25. Interestingly, for a > 0.37 the pairing exists 
only for finite momentum states; i.e. there is a nonzero lower critical momentum at which 
the pairing disappears. The maximal critical values at which the pairing disappears in the 
whole a-P plane are ac = 0.41 and Pc = O.Spp. The main conclusion that can be drawn 
from the discussion above, is that two phase transition take place as the isospin asymmetry 
is increased: first a phase transition from the BCS superfluid state with P = to the LOFF 
superfluid state with P 7^ and, second, a phase transition from the LOFF state to the 
normal (unpaired) state. 

Figure 2 displays the latent heat of phase transition Q = (S'„ — S'^)//? as a function of the 
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isospin asymmetry a and total momentum P. At the boundary of the phase transition from 
superfluid to the normal state in the a-P plane Q = 0, 5"^ = S^] hence the phase transition 
is of the second order (recall that this result holds in the mean-field approximation used in 
determining the entropy). In contrast, Q 7^ 0, for the phase transition from the BCS to 
the LOFF phase and the phase transition is of the first order, except along the line of the 
intersection of Q{pi^ P) surface with the a = = P plane. Note that this line marks the 
region with anomalous negative sign of Q (i.e., in this region the entropy of the superfluid 
state is larger than that of the normal state). This anomaly does not result in a metastable 
state, as the net change of the free energy shown below remains always negative. 

Figure 3 shows the difference of the free energies of the normal and superconducting 
states {SF)pj3, which is normalized to its value in the symmetric and zero-total-momentum 
case ((5Foo)p/3 = -7.35 MeV at p = 0.16 fm'^ and p-^ = 3 MeV. The onset of the LOFF 
phase is seen by examining the constant a slices of the {5F)pp surface. The onset of the 
LOFF phase is signaled by the change of the shape of the these curves: for a > 0.25 the 
minimum of {5F)pp as a function of P shifts from the P = to P 7^ values; i.e., the 
ground state energy corresponds to the state with a total nonzero momentum of the pairs. 
The minimum of the free energy difference, as the maximum of the gap function, is located 
at P ~ 0.3pi? and is not sensitive to the value of a > 0.25. The similarity of the functional 
dependence of the free energy difference and the pairing gap on the parameters a and P is 
not accidental, as {5F)pp is dominated by the pair interaction (condensation) energy given 
by the second term in Eq. (|1^, which scales as pairing gap squared. 

IV. SUMMARY AND OUTLOOK 

In this work we have analyzed the BCS solutions for the neutron-proton pairing in the 
asymmetric nuclear matter when the Cooper pairs are allowed for a nonzero total momen- 
tum. The quasiparticle excitation spectrum is four-fold split compared to the usual BCS 
spectrum of the symmetric, homogeneous matter. The twofold splitting occurs due to the 
finite momentum of the condensate and/or the isospin asymmetry and/or the finite quasi- 
particle lifetime; the simultaneous paring in the isospin single and triplet states leads to a 
further twofold splitting of the spectrum. The gap equation, which was solved numerically 
in the limiting case of vanishing isospin triplet pairing, has nontrivial solutions with finite 
total momentum of the pairs. The corresponding nuclear LOFF phase is the true ground 
state of the system for density asymmetries larger than 0.25. The minimum of the free 
energy corresponds to the total momentum of the condensate P = Q.?)pp independent of the 
value of a. For sufficiently large asymmetries (a > 0.3) the condensate can exist only in the 
state with finite momentum; i.e., apart from a upper critical total momentum for vanishing 
of the condensate, there is a lower one at which condensate sets in. The maximal values of 
the total momentum and asymmetry that the condensate can sustain are P = 0.3^^? and 
a = 0.41. The actual value of a found for the nonrenormalized single particle spectrum 
could be reduced by a factor of 3 if a renormalization is carried in the Brueckner theory of 
nuclear matter. 

Thus, in a definite region of the a-P plane the neutron-proton condensation occurs at 
nonzero momentum of the Cooper pairs, which leads to the formation of a spatially inhomo- 
geneous phase of nuclear matter. This implies a periodic (translationally and rotationally 
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invariant with respect to the basis vectors) spatial structure of the condensate which carries 
an isospin density wave at constant total number of particles. One of the consequences of 
the periodic structure is that the quasiparticle velocities in certain directions could be close 
to zero, which implies a strong anisotropy of the kinetic coefficients of the matter and larger 
heat capacity than in the homogeneous phase. 

The phase transition from the LOFF phase to the normal (nonsuperconducting) phase is 
a transition of the second order. However, the phase transition from the BCS to the LOFF 
phase turns out to be of the first order; i.e., there is a latent heat of transition associated 
with this phase transition. In a certain region of the a-P plane the latent heat has an 
anomalous negative sign. However, this does not affect the stability of the LOFF phase, 
since its energy budget is dominated by the pair-condensation energy. 

What lattice structure prefers the nuclear LOFF phase? The problem of the energetically 
most favorable structure of the LOFF phase has not been solved so far in general. For small 
gaps the integral equation ( |T^ ) is linear and we can seek the solutions in terms of a Fourier 
expansion 

A(r-)=5:A„e^^"^ (20) 

n 

where the lengths of the "basis vectors" are equal. Fulde and Ferrell studied in detail the 
thermodynamics of the LOFF phase with the order parameter containing a single harmonic: 
A(r) = Aqc*"^'' [|l5l. Perhaps on symmetry grounds one can argue that a symmetric ansatz 
A(r) = A(— r), which implies a real gap function, is the case. In the latter case the most 
general form of the harmonic expansion is 

A(r) =2^A„cos(P„r). (21) 

n 

The limiting case of a single harmonic A(r) = 2AoCos(Pr) has been studied by Larkin 
and Ovchinnikov in this case one finds a layered structure. Perhaps, a cubic structure, 
in which case A(r) = 2Ao [cos(Px) + cos{Py) + cos(P2;)], is preferred to the layered one if 
there are no preferred directions in the problem. To conclude, the periodic structure of the 
LOFF phase has been studied only for limited configurations or spatial dimensions so far. 
The determination of the true ground state structure of this phase remains for the future 
work. 
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FIGURES 




FIG. 1. The pairing gap A{pf) in the ^Si-^Di partial wave channel as a function of the isospin 
asymmetry a = {pn — Pp)/p and total momentum P in units of the Fermi momentum. The pairing gap is 
normalized to its value in the symmetric and zero-total-momentum case Aqo = 14.38 MeV at fixed total 
density p = 0.16 fm""^ and temperature (3~^ = 3 MeV. 
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FIG. 2. The latent heat of transition as a function of isospin asymmetry a and total momentum 
P. The remaining parameters are the same as in Fig. 1. 
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FIG. 3. The difference of the free energies of the normal and superconducting state {SF)pf^, 
normalized to its value in the symmetric and zero-total-momentum state ((^-Pbo)p/3 = —7.35 MeV. 
The remaining parameters are the same as in Fig. 1. 
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